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Molecular clouds (MCs) are stellar nurseries, however, formation of stars within MCs 
depends on the ambient physical conditions. MCs, over a free-fall time are exposed to nume- 
rous dynamical phenomena, of which, the interaction with a thin, dense shell of gas is but one. 
Below we present results from self-gravitating, 3-D smoothed particle hydrodynamics (SPH) 
simulations of the problem; seven realisations of the problem have been performed by vary- 
ing the precollision density within the cloud, the nature of the post-collision shock, and the 
spatial resolution in the computational domain. Irrespective of the type of shock, a complex 
network of dense filaments, seeded by numerical noise, readily appears in the shocked cloud. 
Segregation of the dense and rarefied gas phases also manifests itself in a bimodal distribution 
of gas density. We demonstrate that the power-spectrum for rarefied gas is Kolomogorov like, 
while that for the denser gas is considerably steeper. As a corollary to the main problem, we 
also look into the possibly degenerative effect of the SPH artificial viscosity on the impact of 
the incident shell. It is observed that stronger viscosity leads to greater post-shock dissipati- 
on, that strongly decelerates the incident shock-front and promotes formation of contiguous 
structure, albeit on a much longer timescale. We conclude that too much viscosity is likely to 
enhance the proclivity towards gravitational boundedness of structure, leading to unphysical 
fragmentation. On the other hand, insufficient resolution appears to suppress fragmentation. 
Convergence of results is tested at both extremes, first by repeating the test case with more 
than a million particles and then with only half the number of particles in the original test 
case. 
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1 INTRODUCTION clumps, some of which could perhaps spawn stars. This dynamical 

process of forming stars apparently occurs on a short time-scale, 
Expanding shells and triggered star-formation. The stability of „ , . . ..... . , , 

oi about a crossing-time; a similar idea was previously proposed 
the galactic disk and its evolution depends on the availability of ^ lElmegreenl (feoOO) 
molecular gas, that can be converted in to stars. Thus, an important 

prelude to understanding galactic evolution is the determination of _ , , ... „ ,. .,„,,, T , 

Turbulence in the interstellar medium (ISM). Large scale tur- 

factors that could possibly control the rate of star-formation. The , , . , , ■ , • , , , , ,■»„/-. % 

bulence in the galactic disk, and within molecular clouds (MCs), 

problem in turn behoves a study of the stellar feed-back mechanism , ...... , . ... , ■ , 

r , J . appears to be sell-similar and scale-free. Various physical pro- 

and its effect on the efficiency of star-formation. On this latter issue . , , , , , , . ,. - , 

perties such as the length-scale, and velocity dispersion of the 

there has been considerable harrumph leading to antipodal sug- , , . . , , , ■ ■ , 

, I °—— c — 1 1 ~~a observed turbulent motion were related by empirical expressions 

gestions : one such view propounded by | Krumholz & Tan| 120071) first ted b Larson fj^ Turbulent ^ ^ B 

for instance is that, stellar feed-back, and turbulence in general act . , . . , , . „ . .... , , , ■ , 

f inferred via broadening of spectral emission lines, and believed 

to vitiate star-formation as a consequence of which the process , ... , . , , , .. , ~ 

n to be crucial in producing local density enhancements, is often 

may be retarded, resulting in rather inefficient star-formation. , . . , , . . , , „ . .. , .. „ 

J „ r „ thought to be shock-induced. Fluctuations in the density field 

This is reflected m a small value of the star-formation efficiency , , . , . , . , . , , 

produces hierarchical structure, evident from the preponderance 

calculated over a free-fall time (SFEff). This proposition has been - , n . , , . T „, T „ , 

I — LLL fj — L — 1 of dense filaments and clumps in the ISM. Some authors have 



contested by Ballesteros-Paredes & Hartmann ( 2007), who suggest 
that turbulent motion assists rapid assembly of gas into dense 
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previously also employed statistical techniques to study the density 
field, and in this connection, a particularly useful quantity is the 
probability distribution function (PDF) of gas density. The non- 
uniform nature of a turbulent density field for an isothermal gas 
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has been demonstrated to be approximately lognormal with the aid 
of numerical simulations (e.g. Vazquez-Semadeni 1994; Padoan, 
Nordlund & Jones 1997a; Scalo et al. 1998; Vazquez-Semadeni et 
al. 2008). 

In fact, some authors, using the the lognormal density PDF as 
a basis, have derived a general dense core mass function (CMF) 
and the stellar initial mass function (IMF), that is roughly lognor- 
mal with a slope consistent to that of the canonical IMF (e.g. Pa- 
doan 1995, Padoan, Nordlund & Jones 1997b; Padoan & Nordlund 
2002). The latter authors in support of their hypothesis, have also 
sh own that the densi ty distribution in obscure dark-clouds reported 
bv lLadaef a/1dl994l) was consistent with a lognormal density PDF. 
The apparent similarity in the nature of the density PDF with the 
stellar IMF, and/or the CMF, though intriguing, may n ot necessari- 
ly hav e a causal connection as has been pointed out by Scal o et al] 
(1998). Yet others have sought to explain the origin of hierarchical 
structure through the influence of turbulence over fractal MCs (e.g. 
Elmegreen 1997; Elmegreen 2002). Turbulence in the ISM therefo- 
re plays a crucial part in creating density structures, and perchance, 
regulates the rate of star-formation in MCs; shocks supposedly re- 
sponsible for generating turbulent velocity fields could be driven 
either by stellar winds or energetic jets, expanding dense gas shells 
or collision between energetic gas streams (e.g. see reviews by El- 
megreen & Scalo 2004; Mac Low & Klessen 2004, and references 
therein). 

Expanding shells An approximately spherical shell of po- 
werful ionising radiation (the Stromgren sphere), driven by eit- 
her a young star-cluster, a massive O star or a s upernova, sweeps 
up g as in the ISM whence it acquires mass dBallv & Scovillel 
Il98d) . The impact of such a shell, moving supersonically relati- 
ve to the ISM, with a MC could compress the latter, and possib- 
ly enhance star-formation in case the clou d is self-gravitating, or 
trigger star-formation i n a quiescent MC l Lada & Woodenl 1 19761 ; 
iDeherveng et q/.ll2005l ; TZavango et a/.ll2006l) . The MC could even 
be disrupted and sheared apart in case the shock is too strong 
( lElmegreen & Ladj|l977h . Although direct observations of expan- 
ding shells interacting with MCs are relatively scarce, there are 
a few cases where the region of impact has been well studied. 
iRainev et al] d 19871) and other authors cited in that work, using 
emission maps due to shock excited CO molecules in the region 
M17 (the Omega Nebula), have suggested sequential star-formation 
to have been triggered due to t he interaction of an H II shell with 
MCs in the region. Similarly, iKoo & Moon <"1997a.b) have stu- 
died the interaction between a MC and an SNR, W51C using the 
H I 21 cm line and other molecular line emissions. The post- 
collision shock produces filamentary structure in the gas, howe- 
ver, there is no evidence for any star-formation activity. The Vela 
SNR, shown in Fig. 1, is a nother well studied example, see for in- 
stance Nic hols & Slavianl |2004) and references there-in. The ap- 
parently filamentary and other spatially extensive bright spots visi- 
ble in Fig. 1 are supposedly due to MCs engulfed by the SNR. A 
survey of shocked CO emission in the region by Moriguchi et al. 
(2001) has shown that the average size of a MC in the Vela SNR 
region is about 1 pc, having mass of about 60 M Q ; the smallest 
cloud is about 3 M B massive, wit h a size of about 0.3 pc. Recently 
lAzimula. Fich & McCoevI (2009) have reported star-formation ac- 
tivity in clouds S175, believed to be triggered by associated HII 

regions. 

In one of the first numerical studies of the problem, Woo dward! 

dl976h using grid based calculation showed that a supersonic shock 
has a strong shearing interaction with the cloud, that is even- 
tually ruptured by the resulting hydrodynamic instabilities. Simi- 



Figure 1. The Vela SNR photographed by the ROSAT shows several bright 
regions, and an outflow has been detected in one of them. The adjoining 
boxes show respectively, the region where the outflow has been detected, 
and the Chandra image of the outflow. Image adopted from the electronic 
Chandra gallery (http://chandra.harvard.edu/). 



lar results were repor ted by llGein. McKee & Colellal j 1994b and 
Mac Low et al, dl994l) . in case of magnetised clouds . IPongracid 
( 1994) simulated the problem using self-gravitating smoothed par- 
ticle hydrodynamics (SPH) and showed that self-gravity of the post- 
collision shocked gas was crucial for star-formation to commence 
in it; the MC is likely to diffuse away in case gravity remained sub- 
servient to thermal energy, in other words when the post-shock ra- 
diative cooling was inefficient. In the present work we discuss high 
resolution SPH simulations and wish to emphasise on the density 
distribution within the shocked cloud. While the previous work em- 
ploying grid codes has shown that a strong shock flattens the post- 
collision cloud on a rather short timescale, without any other dy- 
namical signature within the cloud, this work suggests, such is not 
the case. We argue that shock induced turbulence leads to structure 
formation within the shocked cloud, some of which shows signs of 
boundedness. Convergence of the results is tested by repeating one 
of the test cases with eight times the number of particles within the 
precollision cloud, and then with half the number of particles. We 
have, however, not tested the extreme case where the MC may be 
flattened by the impact of the incident shell. 

The plan of the paper is as follows. We discuss the initial con- 
ditions in §2 followed by a brief discussion of the numerical sche- 
me §2.1. Various aspects of individual test cases and the general 
evolution of the shocked cloud are discussed respectively in §3, 
and §3.1. The results are presented in §4 with conclusions in §5. 
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Table 1. Summarising physical details of the simulations discussed in the 
paper. The artificial viscosity parameters (a, ft) have also been specified. 



Serial 
number 


Mass of^_ 

individual particle (M Q ) 


Physical parameters^ 
of the cloud 


ISOTHERMAL SHOCK 
{a, ft) =(0.1,0.2) 


1 


m,=1.74xl0- 2 
m ; =1.74xl0~ 5 
m t =3.67xl0~ 5 


M cW =500 Mq, # cH =0.5 pc 
Prid ~ 6 x 10 -20 g cm" 3 , 

r cH =i5K 


2 


m,=1.74xl0" 4 
m ; =1.74xl0- 6 
m t =2.06xl0- 6 


M r u=5 Mm, i?rw=0.5 pc 
p cW ~ 6 X 10~ 22 g cm" 3 , 
r cH =80K 


3 


same as in case 1 


same as in case 1 
0"drf=O.91 km s - ', »,„> ~ 1^ 


4 


same as in case 2 
(ADIABATIC SHOCK) 


same as in case 2 


5 


(a,J3) =(1,2) 
same as in case 1 
(ISOTHERMAL SHOCK) 


same as in case 1 


6 


HIGHER RESOLUTION 
m,=2.20xl0~ 3 
mj=2.53xl0- 6 
m t =1.26xl0 -5 


same as in case 1 


7 


WITHOUT ARTIFICIAL CONDUCTIVITY 
same as in case 1 


same as in case 1 


8 


LOW RESOLUTION 

m,=2.93xl0- 2 
f«j=2.93xl0- 5 
m i =1.09xl0 -4 


same as in case 1 



mi,mj,m^ -belonging, respectively to the cloud, the ICM, and the slab 

^ M dd'^cld'^ ete' v cld*Pcld" ie>s P* x ^ pl,e ty mc niilss o£\Sie cloud, lis radius, temperature, precollision velocity, 

and the average density 

C 

a V( ' r -Virial coefficient of the precollision cloud 



2 INITIAL CONDITIONS 

Figure 2 shows a schematic representation of the cloud-shell sy- 
stem. To compensate for our limited computational resources, we 
consider only a section of the shell, modelled as a thin, dense slab 
having temperature, T shen , and moving at velocity, V s . We admit 
that such a choice for the initial conditions will reduce the post- 
collision momentum delivered to the stationary cloud. The molecu- 
lar cloud (MC) is modelled as a uniform density sphere, confined 
by a diffuse, warm intercloud medium (ICM) [<1 cm 4 ; and a few 
times 10 4 K]; pressure equilibrium is maintained at the cloud-ICM 
interface. In SPH, the ICM is represented by special type of partic- 
les that exert only hydrodynamic force on ordinary gas particles. 
The entire assembly is enclosed in a periodic-box that is used to 
ghost particles, i.e. particles leaving from one face of the box arrive 
from the opposite face. The simulations were terminated once the 
slab traversed the box-width. Six test cases had ~400,000 partic- 
les with ~ 30,000 particles in the cloud, and ~ 80,000 particles in 
the slab; the rest were ICM particles. The average SPH smoothing 
length, h avg , within the cloud is 




A^/,;=30,000; for number of neighbours of an SPH particle, 
N ne n,s-50, hayg ~ 0-03 pc < Aj = (na 2 Q IGp) 111 ~ 0.2 pc, the Jeans 
length for gas at 20 K and average density in the shocked cloud, p ~ 
10~ 19 g cm 4 . We repeated the calculations for one of the test cases, 
case 1 listed in Table 1, with an eightfold increase in the number of 
particles within the cloud alone, with ~240,000 particles in the slab, 
and 730,000 particles in the ICM so that the computational domain 
collectively had 1.2xl0 6 particles. Then for the same of the num- 
ber of neighbours, N„ ei t, s , h mg ~ 1.5 x 10~ 2 pc, and h avg /Aj ~ 13, 
sufficient to resolve the Jeans instability according to the Truelove 
criterion (Truelove et al. 1998). For an arbitrarily chosen mass of 
the cloud, M cU , its radius, R cU , is deliberately chosen much smaller 
than that prescribed by the Larson's relation 

tf drf (pc) = 0.1(pc)(^)° 5 , (1) 

jLarsorJl98ih . The sound-speed, a , inside a self-gravitating cloud 
must be much smaller than its escape velocity, v g , defined as, v 2 = 
1G " c,d . The temperature within the cloud, T c u, is chosen such that 
Oj « vj. The turbulent velocity field i n the MC in case 4 is set up 
in a way similar to that described by M ac Low et al. | Jl997h . The 
velocity perturbations are so set that the initial power spectrum is 
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P{k) oc k° , a = 0; the initial amplitude of the v elocity perturb ations 
is set using the velocity scaling relation due to lLarsor] dl98 lh . The 
physical parameters for each test case have been listed in Table 1 
above. 



2.1 Smoothed particle hydrodynamics (SPH) 

The simulations discussed here have been performed using SPH, 
a Lagrangian, particle based scheme. An SPH particle, contrary to 
intuition, has a finite size characterised by its smoothing length, 
h. Thus, an SPH particle carries physical properties such as mass, 
density, and velocity; SPH particles interact with each other via nu- 
merical viscosity, called artificial viscosity (AV). In all test cases, 
bu t case 5, we have used the att enuated form of viscosity prescribed 
by |Morris~& Monagharj d 1997b . that has been demonstrated to bet- 
ter represent shocks. Under the scheme SPH viscosity parameters 
(a,P) = (0.1,0.2). In addition, we also use the Artificial thermal 
conductivity scheme devised to smo oth out disc ontinuities across 
regions with steep density gradients (|PriceJ|2008). In the fifth case 
we have employed the conventional SPH viscosity, with (a,/3) = 
(1,2), combined with Artificial conductivity. This is intended to de- 
monstrate the detrimental effect of SPH viscosity on shocks, and the 
formation of post-shock structure through the interplay of various 
dynamical instabilities. We used the well tested SPH code, Seren, 
for our purpose here. The robustness of the code has been demon- 
strated by the authors through a suite of tests, that includes model- 
ling hydrodynamic instabilities such as the Kelvin-Helmholtz insta- 
bility (KHI) and the non-linear thin-shell instability (NTSI) (Hub- 
ber et al. 2010 a,b). The code uses a Barnes-Hut tree dBarnes & Hutl 
1986), to find the nearest neighbours of an SPH particle and to cal- 
culate the net force on it. Seren also includes the quadrupole mo- 
ments of distant cells in its calculation of self-gravity. Seren em- 
ploys the modified M4 kernel suggested by Thomas & Couchman 
(1992). 

As in our previous work ( Anathpindika 2009), thermodynamic 
details in the first four cases here have been tackled using a simple 
barotropic equation of state (EOS), defined by Eqn. (2) below, and 
plotted in Fig. 3, while the post-collision shock in case 5 is purely 
adiabatic, defined by the third part of Eqn. (2). (The adiabatic gas 
index, y = 5/3). 



;p < 10- 23 gcirr 3 



- = (k B /m) x 
P 



j ; 10~ 23 g cm" 3 < p < 10~ 22 g cm" 3 

< p < 5 x 10 22 g cm 3 

;5 x 10- 22 gcirT 3 <p 
s: 10~ 18 g cirr 3 

p > 10~ l8 g cirr 3 . 



(2) 



The first two parts of the EOS are designed to maintain the suita- 
ble temperatures in the shell, and the cloud, respectively. The im- 
mediate post-shock temperature within the MC is allowed to rise 
adiabatically using the third part of the EOS after which the gas 
is cooled down to its precollision temperature, T c u, as can be seen 
from the fourth part of the EOS. The gas downstream of the shock, 
after becoming sufficiently dense, is cooled further and maintained 
at temperature, T iso = 10K. The last part of the EOS, that allows 




Figure 2. Schematic representation of the cloud-shell system in the plane 
of collision. A section of the expanding shell, shown as a plane slab, mo- 
ves towards a stationary cloud with a velocity V,. Note that the cloud is 
embedded in a warm intercloud medium. 




Figure 3. The EOS defined by Eqn. (2) above in the logarithmic T—p space. 



the temperature to rise adiabatically when condensation occurs, is 
never broached in these simulations. 




Figure 4. A plot demonstrating the stability of the preshock cloud in isolati- 
on. It shows the probability distribution, p p , of the density of SPH particles 
in a relaxed cube of unit dimensions. The distribution peaks at unity, the 
true density of the box. 
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3 EVOLUTION OF THE SHOCKED CLOUD 
3.1 Stability of the precollision cloud 

Before proceeding to discuss the evolution of the shocked cloud, 
we wish to demonstrate that the observed structure in the shocked 
cloud is indeed physical, produced due to propagation of shock in 
the cloud. The initial conditions for the numerical experiment were 
set up using a relaxed configuration of particles distributed uni- 
formly in a cube of unit length. Particles in the cube were relaxed 
by evolving the box for about 5 sound-crossing times so that its ave- 
rage density is within a few percent of its true value, p trm . Figure 

4 shows the probability distribution for the density of SPH particles 
in the relaxed box, which as expected, peaks near Pt rlte . The requi- 
red uniform density cloud was scooped out of this relaxed box that 
does not show evidence of clumping i.e. spurious density enhance- 
ments, over the relevant dynamical timescale, the crossing time, t c , 
for the cloud, defined by Eqn. (3) above. 



3.2 Summary 

Although the set of test cases discussed in this paper could possibly 
be divided into 3 subgroups as indicated in Table 1 above, the gene- 
ral nature of evolution of the shocked cloud in all cases is rather si- 
milar. The precollision cloud in each case is colder and denser than 
the incident slab, and so, as a direct consequence of the equation of 
continuity, it follows that the post-collision shock would propagate 
within the cloud at a much lower velocity. This, indeed is observed 
in the simulations and the shock-wave propagating within the cloud 
forms a wake, a bow like structure. The impact of the incident slab 
on the cloud surface reflects a shock wave in the ICM, as is evident 
from the rendered density plots of Fig. 5, and the supplementary 
animation. These plots show the cross section of the cloud taken 
through its mid-plane and the density is measured in M pc~ 3 ; ti- 
me on all rendered plots is marked in Million yrs. Gas shearing 
the cloud surface renders it unstable to the Kelvin-Helmholtz (KH) 
instability, while the cloud accelerated by the shock-wave within 
its interiors renders it Rayleigh-Taylor (RT) unstable. The incident 
slab, due to its finite thickness losses momentum after shocking 
the ICM and the cloud, and consequently the pressure at the rear 
end of the cloud is insufficient to trigger its implosion. So, despite 
our observation in the present work being appar ently cont r adicto - 
ry t o that reported previously by, f or instance [Woodward ( 19761) . 
and lKlein. McKee & ColellJ i l 19941) . in fact, there is none; the pre- 
sent work effectively tackles the case of a thin shell impacting a 
cloud which is clear from the argument below in §3.3. We argue 
that shock-induced turbulence generates dense structure on a rather 
short timescale, some of which could become self-gravitating and 
go on to produce prestellar cores. This result holds irrespective of 
the nature of shock, as can be seen from the .comparative plots in 
Fig. 6. Convergence of results is tested by repeating case 1 (§4.2) 
with eight times the number of particles in the preshock cloud and 
then at the other extreme, with approximately half the number of 
particles; see Fig. 12 below. Next, we discuss our results in some 
more detail. 



3.3 Dynamical evolution of the shocked cloud 

Before embarking upon a discussion of various test cases, it will 
be useful to define some timescales which will be summoned in 
further analysis. A slab moving with a velocity, V„ would traverse 
across a cloud of radius, R c u, in a crossing time, t c , defined as 



tc 



2R C , 



(3) 



The cloud crushing time, t cc , the timescale over which a shocked 
cloud might flatten is, 



X U2 2R cU 



(4) 



where y = — — 1000, in all the simulations discussed. 

Of greater importance is the dynamical timescale, t d , of the 
incident shock, the timescale over which pressure behind the in- 
cident shock may vary. When t cc <K td, the shock is practically 
unaffected even after impacting the cloud. On the contrary, when 
td t cc , the dynamical properties of the shock change substanti- 
ally after impacting the cloud. Thus, a qualitative definition of a 
small and a large cloud may be adduced to these preconditions: a 
cloud in which t cc <K t d , may be classified as small, while a large 
cloud is one with t d <k t cc . For illustrative purposes we shall ad- 
opt the connotation of denoting the preshock state variables with 
the subscript '1', and corresponding post-shock variables with '2'. 
The post-shock velocity of the slab can be obtained simply from 
the equation of continuity as, 

V !2 = (— K (5) 

\Pcld2' 

The pressure, p c u 2 , in the shocked cloud is given by the 



pressure-jump condition as, 
Pern 2 7 M 2 - (y - 1) 



(6) 



Pcm (r + 1) 

M is the precollision Mach number (Shore 2007). For a post-shock 
temperature, Ti, the corresponding sound speed, a%, is related to 
the pressure simply by, 



Pddl a iPcU2 



(7) 



Pcm a^fiddi 

The maximum post-shock density for a radiative shock is 
Pcua ~ M 2 p c ui; so that combining this condition with Eqns.(5), 
(6), and (7), we get, 



\pcui >\a 2 l 



Vs2 
V c ld2 



(8) 



also for a strong shock, it follows from Eqn. (6), 

£2£ ~ 2 7 M 2 . (9) 

Pcldl 

Equations (8) and (9) led us to- 

\a 2 l V cU2 

Now, for a shock that has traversed the cloud radius, and re- 
mained fairly stable, R dd ~ t d v cU2 , so that from the equation above 
we have: 



t d 



^2y(^ 2 

V s2 \a 2 l 



(10) 



As a direct consequence of the temperature-jump condition, we ha- 
ve 

/oo\ 2 = (y + I) 2 
\a 2 ) 2jM 2 (j-l)' 

so that Eqn. (10) finally becomes, 

RrM (7+ I)' 



td ~ 



V s2 M 2 (y- 1)' 



(11) 
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Figure 5. A scroll showing the mid-plane, cross-sectional rendered density plots (measured in M Q pc~ 3 ) of the shocked cloud in case 1. The high velocity slab 
impinging on the cloud from the left drives a shock-wave into the latter and injects turbulence within the cloud; shortly thereafter, (f = 0.017 Myrs), turbulence 
produces a network of filaments, much denser than the precollision cloud. Animation available for online version only. 
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Figure 6. From left to right, the panels show a rendered cross-sectional density plot of the shocked cloud on termination of the simulations in cases 1, 3, and 4, 
listed in Table 1 above. It appears that an additional turbulent velocity field assists structure formation, as is evident from the comparatively shorter timescale 
on the plot in the central panel, relative to that on the left-hand panel. The plot in the right hand panel, that of an adiabatic shock, shows less fragmentation 
compared to the earlier two cases. See text for description. 



In the present work, with y = 5/3, and M ~ 25, so 



where r\ ~0.1. 

Cases 1 and 3 (Velocity of the incident shock, V s = 200 km 

The front surface of the slab, on colliding with the stationary cloud 
generates a shock wave in the latter, and after a brief interval of 
time there is a second, weaker shock when the trailing end of the 
incident slab shocks the cloud. This double shock is evident from 
the step in the red curve of the top left-hand panel of Fig. 8 below. 
The initial impact of the slab reflects a rarefaction wave in the ICM 
at the periphery of the cloud, characterised by the leftward shift of 
the green curve relative to the red curve, and a dip in the former near 
the front-edge of the cloud, which after being shocked, has now be- 
come denser. This is also evident from the plot in the central panel 
of Fig. 5, in which the denser edge of the shocked cloud is visible. 
The shock propagating within the cloud however, is progressively 



weakened as evidenced by the progressive drop in pressure within 
the cloud at latter epochs, shown in the top left-hand panel of Fig. 
8. 

The propagating shock generates turbulence with the cloud 
which itself is dissipative, and produces fractal structure which at 
latter epochs, appears more filamentary with a few clumps. The 
left-hand panel of Fig. 6 shows the fragmented interiors of the 
shocked cloud, the velocity vectors overlaid on the plot show the 
underlying chaotic field. Apparently, this structure develops on a 
timescale comparable to the crossing time, t c , defined by Eqn. (3) 



above. The length of the fastest growing mode, /!„„-;, ~ 




and v e ff ~ (a-o + Vcui); v cW2 ~ V sl IM 2 ~ 0.4 km s~', a ~ 0.27 
km s~' at 20 K. Then for a typical post-shock density, p ~ 10~ 18 
g cnT 3 , A, urb ~ 0.1 pc, an order of magnitude smaller than h avg for 
this case. Thus the fragmentation observed in the shocked cloud 
here is likely to be physical. 

Following a shearing interaction between the slab and the ou- 
ter layers of the cloud, its surface is rendered KH unstable. Alt- 
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hough, here we are able to see only a few blobs, broken away from 
the cloud surface; see for instance the last panel in Fig. 6. This 
could perhaps be the case since the shearing surface in these test 
cases is not spatially extensive, following which there is not suf- 
ficient shearing contact. The thin slab surface in contact with the 
main body of the cloud, however, shows a few rolls, but these are 
not particularly well developed as can be seen from the density 
plot, with overlaid contours shown in Fig. 7. Another possible rea- 
son could be the inability of SPH to handle density contrasts within 
fluid layers, which is well documented in contemporary literature. 
While remedies have been suggested to alleviate some of the pro- 
blems related to fluid mixing, their robustness in presence of self- 
gravity is untested. We shall discuss one such exploration in case 7 
below, where the issue of resolving the KH instability in this work 
will be briefly revisited. 

In case 3, we use the same initial conditions as in case 1, but 
the precollision cloud in this case is additionally supported by a 
turbulent velocity field. The amplitude of the turbulent velocity is 
~0.9 km s _I , so that the Virial parameter,av; r , for the precollision 
cloud is approximately unity. This additional velocity field makes 
little difference to the evolution of the shocked cloud. The central 
panel of Fig. 6 shows the fragmented interiors of the cloud, from 
which it is evident that apart from the timescale of fragmentation, 
there is little change otherwise. The cloud in this case appears 
to fragment much faster and develops well-defined structure at 
an epoch somewhat earlier than in case 1. The two cases can be 
compared from the respective plots in the central, and the left-hand 
panel of Fig. 6. It therefore appears that turbulence aids formation 
of structure via fragmentation. 

Cases 2 and 4 (Velocity of the incident shock, V s = 200 km 

s , &y ~2) : The precollision cloud in case 2 is much more ra- 
refied (~ 10~ 23 g cm~ 3 ) as compared to that in case 1, with the 
velocity of the incident slab remaining unchanged. The quantita- 
tive details of the post-collision cloud though, remain identical to 
those of case 1. The initial conditions of this case are also used 
for the fourth case. Having assumed a strongly radiative shock in 
the earlier cases, we adopted an adiabatic EOS in the fourth case, 
rendering the post-collision shock, non-radiative. Consequently the 
shocked gas (~ 1CT 2 ' g cirr 3 ) remains warm at a few hundred Kel- 
vin. We notice that the warm gas, as expected, shows lesser pro- 
pensity towards fragmentation and therefore greater contiguity in 
structure. The reader is referred to the right-hand panel of Fig. 6. 
It is unlikely for such fragments to be bounded and therefore, may 
simply be transitional features. 

A common feature of all these simulations is the relatively 
small value of the dynamical timescale t d , defined by Eqn. (12) 
above. We observe that in all the test cases td <k t cc , so that pressure 
behind the slab changes even before it traverses across the cloud; 
consequently the incident slab is devoid of compressive strength by 
the time it reaches the rear end of the cloud, and so the cloud does 
not flatten. The timescale, t d , calculated using Eqn. (12) above is 
~ 10~ 4 Myr <K t cc ~ 0.2 Myr. The problem, in the form tackled 
here, therefore invariably reduces down to a shock interacting with 
a large cloud or alternatively, a weak shock impacting a cloud. 

We are, however, unable to follow further evolution of the 
filaments and clumps since the periodic box enclosing the system 
has insufficient dimensions, and expansion of the reflected shock 
in the ICM, is impermissible; for mixing between opposite sides 
of the ICM shock must be avoided. The simulations were therefore 
terminated after the formation of dense structure. 



Case 5 : While adopting the normal prescription of the SPH 
viscosity, the initial conditions were kept identical to those in the 
first. Not only does the higher viscosity strongly decelerate the 
incident slab, it also dissipates energy within the slab layers and in 
fact, renders it unstable via the growth of bending modes. Figure 
9 shows a rendered, cross-sectional density plot of the slab-cloud 
system in this case. The wiggles on the slab surface appear similar 
to the features of the thin shell instability, but the initial conditions 
were set up such that the smallest unstable mode was much larger 
than the breadth of the computational box. T his wavenumber 
can b e calculated using Eqn. (14) deduced by iMoriguchi et al\ 
for a thin shocked- shell. So these wiggles are possibly a 
manifestation of excessive dissipation due to stronger viscosity, 
and in fact, the viscous interaction renders the shock within the 
cloud, stationary, as is visible from the central and right-hand 
panels of Fig. 9. 

This is also the case in which the shocked cloud shows the 
least fragmentation, and therefore, the structure within the cloud is 
by and large contiguous. The fastest growing mode, A, llr t, in this 
case will be considerably shortened as the contiguous structure be- 
comes denser which is likely to produce a large number of small 
fragments. This case demonstrates that a numerical study of hydro- 
dynamical instabilities is likely to be corrupted by strong viscosity, 
and potentially misrepresent the situation. We defer the discussi- 
on of cases 6, 7, and 8 to §4.2, following that of the generation of 
vortices in the shocked cloud in §4.1. 



4 DISCUSSION 

We have, in the preceding section demonstrated that the general 
evolution of a shocked cloud is rather insensitive to the nature of the 
post-collision shock. Although the extent of fragmentation in the 
shocked cloud appears to depend on the initial conditions, and the 
post-shock cooling. These re sults are radically different from those 
reported by IWoodward d 19760 . it must, however, be no ted that his 
work and indeed, a latter work bv lMac Low et al\ i 19941) with ma- 
gnetised precollsion clouds neglected self-gravity. Self-gravity has 
been included in the present work, and the evolution of the shocked 
cloud apparently depends on the balance between the gravitational 
force per unit area and the ram-pressure due to turbulent motion. 

4.1 Vorticity in the post-shock medium 

An important point demanding consideration is that of vorticity ge- 
nerated in the post-shock cloud, and the surrounding ICM. It is well 
known that an ideal fluid is free of vorticity, the situation with the 
fluid in these simulations or any other numerical treatment differs in 
two ways : (i) the inevitable presence of numerical viscosity in the 
numerical interpolation renders the fluid non-ideal, and (ii) shocks 
in general are not isentropic. The importance of turbulence in struc- 
ture formation and in general, on the process of star-formation via 
the process of turbulent fragmentation, is well appreciated. It is 
therefore quintessential to look in to the possible sources of tur- 
bulence. 

While an ideal fluid being curl-free, always remains so as a 
consequence of the Kelvin's vorticity theorem, this is hardly true 
in case of real fluids or those of the type considered here. In these 
latter fluids vorticity could be created, destroyed, or may simply 
diffuse from one part to another part of the fluid. To see this, let 
us go through the following set of simple equations: The vorticity, 



8 Anathpindika. S. and Bhatt, H. C. 




0.8 



0.6 




1/1 
c 

-S 

O 



0.7 0. 



Figure 7. A coeval, cross-sectional plot comparing the growth of KH instability (KHI) due to shear between the incident slab and the cloud surface, with, and 
without the Artificial conductivity (AC). The plot in the right-hand panel is the one without AC, and apparently there is not much difference in the two cases. 
We do not observe the rolls associated with KHI; density contours have been overlaid on the rendered image to aid identification of underlying structure. 



w=V X v, for a fluid having velocity, v, the time derivative of which 
is, 

do 

— = V X v. 

dt 

The quantity cj X v = v • (Vv) — 5 Vv 2 , has a vanishing curl, so that 
the above expression may be re-written as, 



do 

— + V X (o) X v) = V X v. 
dt 



(13) 



The equation of fluid motion including the viscosity is, 
£l = _(Ivp + vV 2 v)-V^; 

where v is the kinematic viscosity of the fluid, and (/> g is the gravita- 
tional potential. The remaining symbols have their usual meaning. 
Using this expression in the right-hand side of Eqn. (13), and follo- 
wing a little manipulation we wind up with, 

^ + V X X v) = -[vp X V^ij + vV 2 w|, 

which is the equation describing the time evolution of the vorticity 
of the fluid. The second term on the right-hand side of Eqn. (14) 
above defines the spatial distribution of vorticity in the shocked 
gas. Now using the Laplacian identity, this term can be re- written 

as 



(14) 



[V • (V • a)) - V x V X o\ = V (o 

The gradient of a curl is zero so that the first term on the left-hand 
side vanishes, and we are left with 

-V X V x o) = -V X V 2 v = V 2 cj 

The quantity (V x V 2 v) has been plotted for the shocked gas in the 
cloud, shown in Fig. 10. The shocked gas, in general, shows smaller 
vorticity which suggests, post-shock, vorticity is dissipated. 

An useful quantity to study the temporal evolution of vorticity 
is its circulation, W, defined as 



W 



■X 



(o • dS, 



(15) 



dS being an infinitesimally small area element. This quantity gives 
the flux of vorticity through an area element dS. We have plotted 
the time evolution of the net vorticity in cases 1 and 5. To remind 
our readers, we may recount that starting with identical initial 
conditions and the EOS, case 5 employs the conventional SPH 
viscosity, while the attenuated viscosity, i.e. the Morris-Monaghan 
prescription is employed in case 1. Figure 11 shows a comparative 
plot of the net vorticity for these two cases. It is clear from these 
plots that a larger SPH viscosity generates greater vorticity in 
the post-shock cloud of case 5, although the vorticity in either 
cases decays with time. However, a larger vorticity will lead to 
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Figure 8. Radial distribution of total pressure within the shocked cloud in the first 4 test cases; cases 1 & 2 are plotted in the top left and right panels while 
3 & 4 are plotted in the bottom left and right hand panels. The plots have been made at different epochs of time as the post-collision shock propagates in the 
cloud. The shock following the impact of the incident slab causes steep discontinuity in the pressure distribution at the cloud edges, r ~ 0.5 pc in the plots 
here. Observe that unlike in cases 1 & 3, there is no reflected shock at the cloud surface in cases 2 & 4, where the precollision cloud is much more rarefied. 
See text for description. 




Figure 9. A scroll showing the slab-cloud interaction in case 5. Notice that the incident slab itself, unlike in the first 4 cases, shows signs of instability and 
develops perturbations. Having impacted the cloud, the incident shock soon loses momentum and appears to have become stationary within the cloud; the 
wake in the shocked cloud shown in the right-hand panel is still at its front end. See text above. 



enhanced dissipation of kinetic energy within the system which in 
turn will have grave implications on its dynamical evolution, such 
as on the growth of instabilities and/or fragmentation adduced to 
these instabilities. In general, enhanced viscosity is likely to assist 
fragmentation and result in a number of small, fictitious fragments 
that could corrupt simulations and generate misleading results. 



4.2 Cases 6, 7 and 8 

Testing convergence: Resolution and artificial conductivity 
These two cases are provisioned for additional discussion due to 
two crucial problems that could call into question the tenability of 
the preceding work. 

(i) A pertinent issue in contemporary computational astrophysics is 
the ability of a numerical method to resolve hydrodynamical insta- 
bilities. Both SPH and the adaptive mesh refinement (AMR) codes 
have had considerable difficulty in resolving instabilities such as 
the KHI and the RTI (Springel 2010), constraining the discussion 
to SPH alone, problems such as segregation of regions with steep 
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Figure 10. A plot showing the spatial distribution of vorticity, where the 
quantity (V X V 2 v) has been plotted in logarithmic units over the x—y plane; 
contours for the vorticity have been overlaid on the rendered image. It can 
be seen that shocked gas in general, shows greater dissipation of vorticity 
(t ~ 0.017Myrs). 
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Figure 11. The time evolution of net vorticity in cases 1 (red curve) and 
5 (green curve). Observe that a larger viscosity in case 5, evident from the 
spike at t ~ 1.15/ cr (~ 7.5 X 10 _3 Myrs), produces greater vorticity in the 
shocked cloud. 

density contrast, or clumping of particles have been demonstrated 
(e.g. Agertz et al. 2007; Junk et al. 2010), and partially remedied 
(e.g. Price 2008; Wadsley et al. 2008; Read et al. 2010). While ad- 
ditional dissipative terms to aid fluid mixing, the artificial conduc- 
tivity (AC), prescribed by Price (2008) has been, to some extent, 
successful in alleviating the problem related to segregation of fluid 
layers in purely hydrodynamical simulations, the KHI in this case 
has been demonstrated to grow on a longer timescale (e.g. Hubber 
et al. 2010b); the switch, however, appears ineffective in presence 
of self-gravity as has been discussed above (§3.3, under cases 1 and 
3). Figure 7 is a rendered density plot showing the shearing inter- 
action between the slab and the surface of the cloud. The plot in 
the left and right-hand panels respectively correspond to the cases 
1, and 7, described in Table 1 above; the former employs the AC 
prescription which is turned off in the latter. A visual examination 
alone readily demonstrates the similarity in the underlying struc- 
ture. 

(ii) Another crucial problem is that of particle clumping, the 
so called tensile instability, symptoms of which we had previous- 



ly reported in SPH simulations of dissimilar clouds (see case 3 of 
Anathpindika 2010). Read et al. 2010 have suggested an improved 
kernel along with a prescription for a large number of SPH neigh- 
bours, typically in excess of 400, to improve sampling of the kernel 
core and repel particles with vanishingly small separation. Howe- 
ver, this prescription is likely to increase the smoothing length, at 
the expense of spatial resolution, as suggested by Eqn. (1) above. 
In the present enterprise we have therefore opted to retain the choi- 
ce of 50 neighbours. Also, the physical conditions responsible for 
the tensile instability are unclear. Clumping of SPH particles is a 
rather localised occurrence, that will primarily manifest itself on 
the scale of a smoothing length, the structure reported in the si- 
mulations discussed here, however, is contiguous, spanning a few 
smoothing lengths. We therefore believe that the structure observed 
in the shocked cloud is physical, and corroborate our claim by re- 
peating the simulation in case 1 by adopting a much larger number 
of particles; see Table 1 for physical details. 

The plot in the top panel of Fig. 12 shows a rendered density 
plot of the collision sequence for case 6, the one with the highest 
resolution. The plots are coeval to those for case 1 shown in Fig. 
5 above. It is clear that the overall morphology of the shocked 
cloud, though similar to cases discussed above, shows more 
filamentary structure and clumping in comparison to that in case 1 . 
Enhancement in structure formation is also visible in the density 
PDF plotted for this case in the bottom panel of Fig. 14. As 
expected, this PDF is broadly identical to that for case 1 shown 
in the top left-hand panel of Fig. 13, except in the density range 
~ (10 3 — 10 45 ) cm -3 . The PDF for case 1 in this interval is 
stunted whereas that in the present case has a more pronounced 
peak, which again points to the promiscuous formation of dense 
structure in this case. We adduce this to an increase in spatial 
resolution which in this case is at least by a factor of two as 
described previously in §2. While on the one hand the interface 
between the cloud and the impinging slab is KH unstable, the 
Jeans instability is crucial to structure-formation, and eventually, 
to the formation of prestellar cores in this structure. Resolution 
of this instability is therefore our concern. As reported by several 
authors (e.g. Bate & Burkert 1997; Hubber et al. 2006) in the past, 
insufficient resolution tends to damp the growth of perturbations, 
which in the problem considered here, are purely gravo-thermal in 
nature. The ratio of the Jeans length, Aj, to h mg in the relatively 
lower resolution cases, ft = Ajlh avg ~ 6, is conservative in 
view of the Truelove criterion for resolving the Jeans instability; 
ft :> 8 to satisfy the Truelove criterion. Having increased the 
number of particles within the cloud by a factor of eight (case 6), 
ft = Aj/h m , g ~ 13, so that the Jeans length could possibly be 
resolved. 

Case 8 - Lowest resolution 

The number of particles within the cloud for this case are re- 
duced by a factor of two, so that the average smoothing length, 
h avg , defined in §2 above is ~ 1.25 times that in case 1 so that 
ft = 7— — ~ 4, and therefore, obviously violates the Truelove 
criterion to resolve Jeans instability. The shocked cloud lacks 
sufficient spatial resolution and consequently shows hardly any 
signs of fragmentation which is not surprising at all. The bottom 
panel in fig. 12 shows a cross-sectional plot of the shocked cloud 
in this case. 

These occurrences at the two extreme choices of N c u, adum- 
brated respectively, by cases 6, and 8, suggest that poor spatial reso- 
lution tends to suppress fragmentation so that the number of partic- 
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les in the computational domain should be increased as much pos- 
sible. The point has been emphasised by numerous authors in the 
past (e.g. Agertz et al. 2007; Commercon et al. 2008, and several 
references therein). An associated problem is of the stability of the 
preshock cloud against perturbations. The cross-sectional plots of 
the shocked cloud shown in Figs. 5, 9, and the top panel of 12, for 
instance, show that density structure appears at the rear end of the 
cloud even before that part is engulfed by the impinging shock. We 
therefore conclude that instabilities observed in the shocked cloud 
are triggered purely by white noise, and it appears, one may have 
to use at least a billion particles to possibly ensure stability. We dis- 
cuss the stability of the density field in Appendix A below where a 
conservative lower-limit on N c u has also been derived. 



4.3 Density PDFs and the core mass function 

The PDF of a variable provides information about its distribution in 
sub-intervals spanning the whole range. The PDF in each of the five 
cases discussed above are distinctly bimodal which indicates that 
part of the underlying gas is distributed in clumps and filaments, 
denser than the rest of it; see figs. 13 and 14. The remainder of the 
gas is unbound and therefore rarefied. The volume filling factor, /, 
of the dense structure within the shocked cloud can be roughly esti- 
mated by taking the ratio of the observed peak in the density PDF 
with the density of the precollision cloud. This yields / ~10 %, in 
other words 90 % of the shocked gas ends up in rarefi ed holes, or 
voids, which is in agreement with the value reported by Elmegreen 
dl997h for fractal clouds. This observed segregation of the dense 
and the rarefied medium is an important characteristic of the bi- 
stable model of the interstellar medium. While this is true for the 
denser precollision clouds in cases 1,3, and 6, the rarefied clouds in 
cases 2, and 4, after being shocked, show some what contradictory 
behaviour. The PDF in case 2 is peaked at a density close to the 
density of the precollision cloud, which suggests, the isothermal 
shock has had little effect on the distribution of precollision densi- 
ty within the cloud. On the other hand, the adiabatic shock in case 
4 leads to a more significant density enhancement in the shocked 
cloud. The shocked gas in this latter case is warm, at roughly a few 
hundred Kelvin which is more than an order of magnitude larger 
than the gas temperature in case 2. The pressure-distribution wi- 
thin the shocked cloud plotted in the bottom right-hand panel Fig. 
8 shows that the cloud in this case has shrunk by about 10% which 
is also responsible for density enhancement within the cloud that 
results in the PDF plotted in the bottom right-hand panel of Fig. 
4. While the incident shock raises the pressure within the cloud, it 
drops off rapidly near edges. 

Fragmentation in warm gas is likely to be delayed, leading to 
greater contiguity of structure. We note that, the PDF in these test 
cases, at the high density end, shows an approximate lognormal 
nature. Our work therefore reinforces the importance of interstellar 
shocks in the process of stellar birth as has been pointed out by nu- 
merous workers in the past (e.g.Vazquez-Semadeni (1994); Padoan 
& Nordlund (2002); Vazquez-Semadeni et al. (2008); Hennebelle 
& Chabrier (2008) ). We need not be unduly concerned about the 
tail of the PDFs extending into spectacularly low densities, as it is 
likely to be part of the warm ICM. At this point it would therefore 
be useful to briefly consider the spectrum for the two phases. The 
energy density of the turbulent field in Fourier space, E(k), and that 
in the real space is related as, 

f E(k)dk = 1 V v 2 , 
Jo 2 i 



k being the wavenumber. The summation on the right-hand side 
above runs over all SPH particles. The smallest spatial scale on 
which energy could possibly be dissipated is the SPH smoothing 
length, hi, so that ~ The spectrum so calculated is plot- 

ted in Fig. 15 which also shows a segregation of the two gas pha- 
ses. While the rarefied gas produces a Kolmogorov like spectrum, 
oc k~ 17 , that for the dense gas is much steeper which suggests that 
the Kolmogorov approximation of inviscid fluids breaks down in 
high density regions. 

We close this discussion on PDFs by briefly considering the 
effect of higher viscosity in case 5, on the distribution of post-shock 
density. As in the previous cases, this PDF is also bimodal, but pea- 
ked at the lower end of the density distribution. This again suggests 
that there is greater contiguity of structure, that may fragment only 
after becoming sufficiently massive. This PDF, as it evolves in time 
will gradually move towards higher densities, however, an import- 
ant constraint here, and particularly, from the perspective of star- 
formation, is the timescale over which this shift may be effected. 
The PDF shown in the top panel of Fig. 14 is plotted at a similar 
epoch as in the previous cases, and it appears that a larger viscosi- 
ty tends to slow down the process of turbulent fragmentation. This 
leads us to two possible scenarios : (a) Turbulence within gas tends 
to produce substructure on a rather short timescale, and (b) quie- 
scent gas on the other hand, tends to agglomerate more gas and 
become sufficiently massive before it could possibly become Jeans 
unstable. The post-collision shock in case 5 becomes more or less 
stationary due to higher viscous dissipation, and so gas within the 
cloud remains by and large quiescent, qualifying for category (b) 
above. The proposition under (b) is less likely to explain the forma- 
tion of stars for it will tend to delay the onset of the star-formation, 
returning a timescale grossly inconsistent with that commonly re- 
ported by surveys of young star-forming regions. 



5 CONCLUSIONS 

Starting with a relatively stable MC, we have demonstrated that tur- 
bulence injected due to the impact of a thin shell produces structure 
within the post-collision cloud. And crucially, this fragmentation 
occurs on a rather short timescale, comparable to the local crossing 
time, t c , defined by Eqn. (3) above. Interestingly, the filling fac- 
tor, /, for dense structure is consistent with that reported for fractal 
clouds in the interstellar medium, which reinforces the idea that 
interstellar shocks could be vital in producing the observed den- 
se structure in the galactic disk; prestellar cores, under appropriate 
physical conditions, could condense in this structure. The degree of 
fragmentation, however, depends on thermodynamic details of the 
process, and in general warm, and/or rarefied gas, as expected, is 
observed to be quiescent. The results of our work are particularly 
interesting in view of the recently reported sub-mm observations 
of the Gould belt with the Herschel space observatory that reveals 
a number of prestellar cores in several dense, gas filaments (e.g. 
Andre et al. 2010). We have tested the convergence of our simula- 
tions by repeating the test with more than a million particles, whe- 
re filaments in the shocked cloud are even more promiscuous; see 
top panel of Fig. 12; insufficient resolution, however, suppresses 
structure formation as can be seen in the bottom panel of Fig. 12. 
Starting with a relatively stable molecular cloud, it has been de- 
monstrated that runaway of growth of perturbations in the density 
field of the cloud lead to formation of small clumps, and more con- 
tiguous filaments. Any density fluctuations arise purely from white 
noise. 
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Figure 12. Top panel A montage showing cross-sectional rendered images of the shocked cloud from the best resolved simulation (N tot = 1.2 X 10 6 ), case 
6 in Table 1. The plot is coeval relative to that shown in Fig. 5 above, for case 1. See text for further description. Bottom panel A montage showing a density 
cross-section of the shocked cloud in case 8. This is the case with lowest resolution and shows hardly any evidence of fragmentation (N ta t = 80, 000). 
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Figure 13. The density probability distributions for the first 4 cases. The bimodal nature of PDFs shows that the propagation of a shock leads to segregation 
of dense and rarefied gas phase. See text for description. 



Turbulence in the physical world is complex, and unlikely to 
be adequately represented by a single power-spectrum, such as the 
Kolmogorov spectrum, and more crucially, the Kolmogorov analy- 
sis was prescribed for inviscid fluids whereas real fluids are viscous. 
Models such as the one tested here provides naturally a more reali- 



stic turbulent field with vortices which permits a study of the forma- 
tion of dense structure from stable initial conditions. Vortices appe- 
ar to play a key role in viscous dissipation of energy, evident from 
the plots in Figs. 1 1 and 12. The montage in the top panel of Fig. 12 
in fact demonstrates that substantial energy dissipation, manifested 
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Figure 14. Top panel : The density PDF for case 5. This PDF, as for the 
first 4 cases, is also bimodal but peaks at a much lower density although the 
epoch of this plot is the same for all cases. Bottom panel : The density PDF 
for case 6, which like other cases also shows a bimodal distribution, but is 
more pronounced for intermediate densities in comparison to that for case 
1 plotted in Fig. 13 which demonstrates preferential distribution of matter 
in contiguous structure. 



by reduced vorticity, precludes structure formation, the bounded- 
ness of which, of course, depends on the extent of dissipation. We 
might thus, be in a better position to answer crucial questions rela- 
ted to the rate and efficiency of star-formation, all of which depend 
directly on the ambient conditions in the star-forming regions. The 
post-collision shock readily generates a power-law distribution of 
interstellar gas as has also been demonstrated by numerous authors 
in the past, while in the cases discussed above, segregation bet- 
ween the dense and rarefied phases of the gas is also evident from 
the PDFs plotted in Figs. 13 and lower panel of Fig. 14. This is fur- 
ther corroborated by the energy power spectrum for the two phases 
plotted in Fig. 15. The spectrum for the diffuse gas is Kolmogorov 
like, while that for the dense gas is much steeper which agrees with 
our observation in an earlier related work (Anathpindika 2010). 

The model tested here, based on the well-known concept of 
turbulent fragmentation, suggests that the process leading to the as- 
sembly of a prestellar core is rapid, and may not last longer than 
a crossing time. The filaments and clumps produced in the present 
work could have been allowed to evolve further. But for the limited 
computational resources at our disposal, any such ideas will have 
to be deferred for future exploration. As a divestment to the central 
problem, we have also attempted to study the influence of nume- 
rical viscosity on the post-shock dynamics. We conclude that too 
much viscosity, and in fact the standard SPH viscosity, is likely to 
produce misleading results by inducing spurious fragmentation. 
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Figure 15. The power spectrum for the two gas phases. Respective power- 
laws have also been plotted for comparison; the power-law for diffuse gas 
shown by a green line, oc k~ 1J , is in agreement with the Kolmogorov spec- 
trum, while that for dense gas, is much steeper shown by a purple line, 
oc k~ 4s . 
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W(r, h) = 



21 



256tt& 3 



(2 - s) 2 (4 - s 2 ); s < 2 



(A2) 



Thomas & Couchman (1992). The time-derivative of the SPH den- 
sity, p, is simply 



dpi ^ dWij(s,h) 
dt ^ dt 



and 



(A3) 



dWjj(s,h) _ dW^fc _ V fj dWjj 
dt ds dt hi ds ' 

which on substituting in Eqn. (A-3) we get 



dpi 



It ~ Ti^ Vii ~dT 



ds 



(A4) 



which is simply the equation of continuity, and v ;j - is the signal ve- 
locity between the particle pair Our interest lies in exploring 
the possibility of structure formation due to a small perturbation in 
the SPH density field so we adopt the standard procedure where the 
perturbation is assumed to be oscillatory with angular frequency o, 



Pt=py io ". 



(A5) 



p° = 21m is the true density of the SPH particle (s = in Eqn. 

' 256nh . 

i 

(A-2)). Making necessary substitutions, Eqn. (A-4) becomes 

21 

—iutp (cos(<uf) - isin(<w0) = m ) v,,(3s 2 - 4), 

64nh 4 . y 

and comparing real parts on either sides we wind up with, 
-to sui(gj*) = 4v;;(3s 2 - 4) 

which for an infinitesimally small perturbation becomes 

-a> 2 t ~ 4v ; ,(3s 2 - 4) 

or 

r v iy (4-3s 2 )ii/2 

«j~2[ ] . (A6) 

This expression leads us to the following three cases : 

Case 1 up- > 0; there exists a real root so that a density per- 
turbation defined by Eqn. (A5) will remain oscillatory. 

Case 2 o) 2 < 0; the root is imaginary which suggests that 
any perturbation, according to Eqn. (A5), will grow exponentially 
in time. 



APPENDIX A: STABILITY ANALYSIS OF THE SPH 
DENSITY FIELD 

Let us begin with the equation for SPH density of the i' h particle, 
piiri), having position vector r, 

Pi = pto) = Wij(s, hi), (Al) 

J 

where s = |r,- — r 7 -|//i,-, W(s,hi) is the smoothing kernel and the 
summation is over number of all SPH neighbours, N ne i bs . Indivi- 
dual particles are assumed to have identical masses, m. The kernel 
employed in the simulations discussed above is 



Case 3 io 2 = 0; the condition when the density field is like- 
ly to remain steady. 

Equation (A6) may be used to derive a lower limit on the 
maximum number of particles, N c u, in a cloud. The condition 
to avoid growth of density perturbations defined by case 3 above 
implies s = 2h 1/2 , or h = (0.5r iy ) 2/3 but r y oc (Ncu)- 1 ' 3 so 
that h oc (N cU )~ 2/9 . Even by the most conservative estimates, 
N c u ~ 10' before the average smoothing length is reduced by 
about 29% to satisfy the equality. Evidently, any possible gain in 
stability will exact enormous computational expense. Thus, for all 
practical purposes the initial set of SPH particles, as qualified by 
case 1 above, is likely to be susceptible to numerical instability. 
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